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Abstract 

In previous work we have established that a particular quasi-classical theoretical view- 
point about the nature of gravitation in galaxy discs provides extremely high concordance 
(comparable with MOND) with the observed dynamics in optical discs. 

In this paper, we consider the possibility that the optical disc and the flat radio disc 
represent distinct dynamical regions which are inherently described by distinct phases of 
the same over-arching theory of which the quasi-classical theory is but one phase. This is 
a defining moment for the theory under discussion since the only natural generalization 
available is that from the existing quasi-classical (one-clock) model into a quasi-relativistic 
(two-clock) model. The quasi-relativistic disc turns out to be a stationary hyperbolic disc 
within which information propagates along characteristics and within which the rotation 
curve is necessarily exactly flat. 

The transition between the two components of the disc is managed by physically de- 
termined jump conditions which are represented by a set of non-linear algebraic equa- 
tions. These equations predict the existence of a forbidden region within the parameter 
space which coincides exactly with that part of the parameter space in which there is an 
already observed highly significant underpopulation of galaxy discs. 

1 Introduction 

1.1 A globally inertial material vacuum 

This is the third paper in a series which has as its purpose a reformulation of gravitation 
theory on the basis of what might be called an 'ultra-strong Machian worldview'. The initial 
paper (Roscoe Q, astro-ph/0107397) asked the question: 

7s it possible to conceive a globally inertial space & time with is irreducibly associated with 
a non-trivial global mass distribution and, if so, what are the properties of this distribu- 
tion? 

This question received a positive answer subject to the conditions that: 

• the global mass distribution concerned is fractal, D = 2; 

• the constituent particles of this distribution are photon-like in the sense that they have 
a velocity property, v, associated with them satisfying |v| = const; 
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• the velocity property is not a true velocity but is, rather, a conversion factor between 
times scales and length scales exactly like Bondi's interpretation of light velocity, c. 

For obvious reasons, the foregoing was considered to be a description of a rudimentary ma- 
terial vacuum. 

1.2 Gravitation as perturbations in the material vacuum 

The second paper (Roscoe H, "astro-ph/0306228) was concerned with the question of how 
gravitational processes could arise in this material vacuum. It began by showing how point- 
wise perturbations of the D = 2 distribution gave rise to classical Newtonian gravitational 
processes for test particles placed in the region of the perturbation and went on to derive the 
equations of motion for an idealized spiral galaxy. This theory was then applied to model the 
rotation curves for a small sample of low surface brightness (LSB) galaxies giving results that 
were directly comparable in quality to those obtained from the MOND algorithm - and far 
superior to anything achieved by the multi-parameter dark matter models. 

1.3 Flat rotation curves and an unsought-for bonus 

However, whilst the theory of the second paper, which we shall refer to as the quasi-classical 
(one-clock) theory, works extremely well within the optical disc, it is not able to generate the 
flat extensions to rotation curves usually associated with the radio disc. The purpose of this 
third paper is to show how a natural generalization of the quasi-classical one-clock theory 
into a quasi-relativistic^ two-clock theory solves this problem in a direct and simple way. In 
essence, we find that the rotation curve of the quasi-relativistic disc is necessarily everywhere 
flat. Thus, in the presented worldview, galaxy discs exist in two phases - one, the optical disc, 
which is described by the quasi-classical theory and the other, the radio disc, described by the 
quasi-relativistic theory. The transition from the quasi-classical disc to the quasi-relativistic 
disc, which can be considered analogous to the transition between a subsonic and a super- 
sonic flow in a compressible fluid, is managed by a set of physically derived jump conditions 
which are represented by a set of non-linear algebraic equations. 

As an unsought-for bonus, it transpires that the jump-condition equations become degener- 
ate in exactly that part of the disc galaxy parameter space which we already know to be very 
significantly underpopulated in such objects in the four large samples that we have analyzed 
(Roscoe [30]). In other words, the theory predicts the difficulty of forming galaxy discs in 
exactly that part of the parameter space which is already known to be significantly under- 
populated by disc galaxies. 

2 Why bother? Some difficulties faced by the standard model 

Although not yet widely accepted as problematical, it is a fact that there exist various well- 
founded observations and difficulties which, if taken at their face value, would cast reason- 
able doubt about the long-term viability of the standard cosmological model. These fall into 
two broad categories which we can classify as redshift phenomena and as non-redshift phenom- 
ena. The point can be illustrated by reference to just two examples which can be considered 
to serve as 'straws in the wind'. 

^By quasi-relativistic, we mean a formulation in which the invariant line element has the general structure 
ds^ = f{x, y, z, t, dx, dy, dz, dt) but within which ds is not defined as an element of proper time. This latter 
quantity arises in an entirely natural, but quite different, way. 
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For more than two decades, Tifft (see |I33| for an early reference) has been making claims about 
the existence of what he terms as redshift periodicities. The claims have become increasingly 
elaborate over the years but, at their most simple, they relate to the supposed existence of 
a periodicity at around 72km/sec between galaxies in clusters and one at around 36km/sec 
between our own galaxy and galaxies in the local group, out to about 2000km/ sec. Unfor- 
tunately, Tifft was never able to support his hypotheses with statistical analyzes rigorous 
enough to match the extreme nature of the claims. However, sometime in the early 1990s 
Napier brought scepticism and a long experience in dealing with noisy periodic data in the 
geophysical and cometary record (variously [?,lll|5l|6|) to the problem . There followed a short 
series of little-known papers by Napier & Guthrie ( l23ll2^ITHl '). These papers were remark- 
able on two counts: firstly, for their statistical rigour and secondly for their conclusion - that 
the Tifft hypotheses relating to the 72km/ sec and the 36km/ sec periodicities were supported 
on the (non-Tifft) samples used by these authors at the level of virtual statistical certainty (see 
specifically 114J). The basic technical problems associated with these analyzes revolve around 
the difficulties involved in finding line-profiles that are symmetric enough to allow HI deter- 
minations of galaxy redshifts to within better than 5km/ sec - it can be shown that anything 
much worse than this makes it impossible to reliably detect any periodicities at the 36km/ sec 
level which might exist. 

So far as non-redshift phenomena are concerned, we need only cite the problems associated with 
modelling rotation curves in disc galaxies in general - and in low surface brightness galaxies 
in particular - using the multi-parameter dark matter models, and then compare this with 
the simplicity and effectiveness with which the one-parameter MOND recipe (see Milgrom 
II20I I2TI l22| for early predictions) succeeds in this task. Of course, MOND is not the answer 
for it is not even a theory - but the details of its successes over a wide variety of applications 
(see variously McGaugh, de Blok and Sanders IITHl[T9ll32l ) should be raising serious questions 
about the standard model's failures in this respect - notwithstanding well-known problems 
experienced with MOND when modelling galaxy clusters. 

2.1 An historic perspective on the proposed alternative view 

Given that the future of the standard model is by no means secure, then there is room to con- 
sider alternatives, such as the one being considered here. This alternative view was originally 
driven by what could be termed as an 'ultra-strong Machian worldview' and this remains as 
the core philosophy. But, interestingly, the theory which has emerged - which employs all 
the machinery more usually associated with curved spacetime theories - falls directly into the 
category of LeSage-type theories. See the introduction of Roscoe |i2J (astro-ph/0306228| for a 
brief discussion of this. 

LeSage 115J was interested in the causal mechanisms which he believed must underly Newto- 
nian gravitation theory and he proposed the existence of an isotropic and homogeneous sub- 
stratum (in modern terms, a material vacuum) consisting of microscopic particles in states 
of agitated motion. These particles were envisaged as colliding with, and transferring mo- 
mentum to, ordinary material bodies. For one body alone in this agitated substratum, the 
overall effect of the transferred momentum was seen to be neutral. However, when two such 
ordinary bodies were near each other, then they provided mutual 'shade' to each other from 
the bombardment of the substratum particles thereby creating an asymmetric imbalance in 
the rate of momentum transfer from the substratum to the two bodies - the net effect being 
that the two bodies experienced net forces which would appear as mutual attraction. When 
the details of schemes such as this are worked through, it is a simple matter to recover a 
quantitatively correct 'Newtonian gravitational force'. 
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The LeSage theory, and many others of a similar kind that have appeared over the years (see 
for example, Edwards [13J for a good review) envisage the substratum particles in a deeply 
classical way - that is, as essentially microscopic versions of macroscopic classical bodies ex- 
isting within a pre-determined inertial space & time. By contrast, the 'substratum particles' 
associated with the theory of the inertial frame developed by Roscoe [IJ (astro-ph/0107397) 
are the very antithesis of the classical momentum-bearing particle: specifically, they do not 
exist within space & time - rather, what we call space & time can be considered as a metaphor 
for the relationships which exist between these vacuum particles - in effect, they define what 
we call space & time. Thus, according to this theory, a globally inertial space & time is irre- 
ducibly associated with a material vacuum within which the vacuum particles are fractally 
distributed with D = 2. Conventional gravitational processes arise in a given locality when 
this distribution is perturbed in that locality. Consequently, an intrinsic assumption of our 
approach is that conventional bodies act as sinks for the vacuum particles and thereby act as 
local perturbers of the D = 2 equilibrium distribution. 

3 Review of the middle-disc data analysis 
3.1 Disc galaxy rotation curves 

When one studies rotation curves 'in the whole' one is typically confronted by various com- 
plexities which would appear to defy any chance of a simple description. These vary from 
non-monotonic behaviour (typically) near the central regions, through generally monotoni- 
cally rising sections through, finally, abrupt transitions into almost perfect flatness somewhere 
in the boundary between the optical and the radio parts of the disc. Astronomers have, tra- 
ditionally, fitted very complicated functional forms to rotation curves in attempts to provide 
global classifications of their shapes, (cf Courteau [7J, Rix et al Il26t ). In this paper, we de- 
scribe the second step of a quite different strategy whereby we consider the possibility that 
galaxy discs might reasonably be partitioned into three distinct dynamical regions, each to be 
considered separately within the context of the same general theory and merged into a whole 
at their joining boundaries by appropriate jump conditions: 

• an innermost region dominated by the central bulge where the dynamics frequently 
appear messy and complicated; 

• a mid-region comprising most of the optical disc where the dynamics appear to be rea- 
sonably regular; 

• an outer region comprising the radio disc where the dynamics gives rise to the very 
peculiar constancy of rotational velocity. 

The phenomenology of the optical mid-regions has already been comprehensively analyzed 
in Roscoe [27, 30 1 whilst the related theory has been developed in Roscoe [1 , 2j (astro-ph/ 0107397^ 
[astro-ph/ 0306228 ). Correspondingly, our main task here is to show how the flat radio disc can 
be understood from the same general theoretical perspective and finally merged seamlessly 
with the optical mid-region. The process by which this is done also makes clear how the over- 
all task, that of including the three dynamical regions within a unified theoretical perspective, 
can finally be completed. 
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3.2 General comments 

The mid-1990's was an optimal time for our partitioning approach since Persic & Salucci 
||25| had just made available a sample of 900 folded optical rotation curves (ORCs) which 
originated from the Southern Sky survey of Mathewson et al 1.16] . This sample was notable 
for three things: 

• being purely optical it contained no radio disc component; 

• it was sufficiently large that it could support a rigorous statistical analysis; 

• the individual ORCs were very accurately folded - this quality being driven by Persic & 
Salucci's intention of using them for studying the interior dynamics of individual discs. 

The availability of this large high-quality sample meant that there was a reasonable chance of 
gaining a genuine insight into the generic behaviour of dynamics in the middle-disc regions 
of galaxies. Furthermore, the quality of Persic & Salucci's folding process turned out to be 
extremely important for it added crucial precision to the statistical study we had in mind. 

3.3 Summary of previous results 

The most simple solution provided by the quasi-classical theory when applied to modelling 
an idealized disc stated that rotation velocities followed power laws V = AR"" where {A, a) 
are parameters that define individual idealized discs. This solution turns out to be merely 
the most simple in a class of much more complex solutions, but it defined the basic question 
which has driven most of this author's data analyzes: specifically, to what extent can rotation 
velocities in real galaxy discs be said to be described by simple power laws? This question turned out 
to be an immensely productive one. The results arising from pursuing it - both directly and 
indirectly - are listed below. 

• Over large samples of late-type galaxies {Sb..Sd), and to an extremely high statistical 
precision, rotation velocities in the outer regions of optical discs can be represented by 
power-laws of the type 



where the scaling constants, {Vo,Ro) are very strong functions of galaxy luminosity 
properties (Roscoe 1127113(1 ^: 

• The power-law structure in the outer regions of optical discs gives automatic rise to 
the classical Tully-Fisher relations which are properly calibrated for J-band luminosity 
data and -R-band luminosity data according to whether the sample being analyzed uses 
/-band or i?-band photometry (Roscoe |301); 

• The pursuit of the power-law question led directly to the discovery of the discrete dy- 
namical states phenomenology in disc galaxies according to which the distribution of 
the parameter A in the law V = AR"' is multi-peaked with four strong peaks satisfying 
InA ;^3.90, 4.17, 4.72 and 5.11 being identified. This translates into the statement that 
rotational velocities at Ikpc tend to cluster around the values of 49km/ sec, 65km/ sec, 
112km/sec and 166km/sec (Roscoe [28. 30J); 
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• The statistical success of the simple power-law formulation for outer-disc optical ro- 
tation curves prompted a much deeper study of the underlying quasi-classical theory 
Subsequently, we were able to show that it could be applied to model the dynamics 
of low surface brightness galaxies (LSBs) giving results that hitherto have only been 
matched by MOND (Roscoe astro-ph/0306228). 

It is the general success of the quasi-classical model, as outlined above, which has encouraged 
the development of the work which forms the main body of this paper. 

3.4 Summary of new results 

The current paper provides three additional significant results which can be stated as: 

• A fresh, but simple, insight into the quasi-classical model provides a qualitative un- 
derstanding of why the power-law model V = AR"" fits the generality of middle-disc 
dynamics as well as it does - even though it is a simple matter to find examples of indi- 
vidual discs which are in clear departure from any power-law description. 

• the rotation curve flatness problem is resolved within the theory: specifically, the flat 
exterior sections of rotation curves arise as the only solutions of the two-clock quasi- 
relativistic generalization of the one-clock quasi-classical theory; 

• the jump conditions between the quasi-classical and the quasi-relativistic parts of the 
disc form a set of non-linear algebraic equations which become degenerate in that part 
of the parameter space which corresponds exactly to the deep valley between the lnA = 
4.17 and the lnA = 4.72 peaks of the In A frequency diagram. That is, the jump condi- 
tions predict the existence of a forbidden zone in the parameter space which appears to 
exist on the data. 



4 Review of the quasi-classical disc theory 

A brief review of the theory underlying the quasi-classical middle-disc model (Roscoe HI 13 
|astro-ph / 0107397, astro-ph/ 0306228 ) provides a useful precursor for the main analysis of this 
paper and provides the first new result. 

4.1 Physical time and mathematical closure 

The basic relation of the quasi-classical theory is given by 

r)2 A// I^M 

where the indices take values 1..3 and A4 provides a scalar representation of the mass distri- 
bution concerned - but see Roscoe [21 (astro-ph/0306228) for a complete discussion. In the 
present case, we are concerned with modelling disc galaxies and here Ai (r) has the straight- 
forward interpretation as the amount of mass contained within a radius r in a disc of unit 
thickness. Furthermore, for disc galaxies, the assumption of cylindrical symmetry is appro- 
priate implying F^^ = 0. Using cylindrical polars, we then find that the line element can be 
written as 

ds = Qijdx dx^ = -^—r-dr^ + — -— r^dO''. 

\r or ' 
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Introducing the quantities R = r/r^ and T = t/to where (ro, to) are scaling constants, we get 
the dimensionless form 

ds" = gijdx'dx^ = ^^dR" + f = AdR' + B R^dO^. (2) 

The equations of motion can now be written down in the usual way using the variational 
principle 

J = y ^/g~¥i^dT. 

However, because I is homogeneous degree zero in the temporal ordering parameter, r, then 
the resulting equations of motion are invariant wrt arbitrary transformations of this parame- 
ter. This has two related consequences: firstly, it implies that r does not, in general, represent 
physical time and, secondly, the resulting equations of motion form an incomplete set. In fact, 
because the analysis is two-dimensional, the variational principal only gives one independent 
equation of motion, namely 



A = R'OB = (^AR^ + BR^e^y 



which corresponds to the conservation of angular momentum in purely classical theory. The 
situation is, of course, identical for General Relativity - there, the equations of motion arising 
from the variational principle are completed by introducing the concept of particle proper 
time. In the case under discussion (see Roscoe 12J or astro-ph/ 0306228 for a complete dis- 
cussion) this was not an option, and a different solution had to be sought: Certain symmetry 
arguments led to the conclusions that, within the context of the application of the theory to 
disc galaxies, 

• the forces acting within the disc are not necessarily central; 

• the ratio of the magnitudes of the transverse forces to the radial forces along any radius 
within the disc is constant. 

This led to the closure of the system in the following way: the line element is given by 

ds^ = AdR^ + BR^de^ (3) 

where A and B are defined in (|2||. Since the closure information is derived from the ratio of 
force components, it is necessary to define what we mean hy force in a disc with line element 
given by (|3|l. To this end, we factorize this latter relation as 



ds'^ = ds-ds, ds = SiVAdR'R + S2VB Rdee (4) 

where 5*1 = ±1, 6*2 = ±1 and where, in the standard notation, R = (cos sin 6^)^ and 
9 = (— sin6', cos^^)-*^. Hence, introducing an arbitrary temporal ordering parameter, r, we 
can define a generalized velocity as 

V = SiVarr + S2VB Roe 

from which a generalized acceleration vector can be derived as 

V= jsi^ (yiR^ -S2VBRe^'^ R+ 1^2^ {yBRO^ + SiyfARO^ 6. 



However, we also require that the ratio of the magnitude of the transverse force component 
to the radial force component is to be constant, so that we get the additional equation 

jsi^ (yiR^ -S2VbR9^'^ =k!^S2^ (yBR9^ +SiVaR9^ (5) 

for some constant K. Since this equation is not invariant under arbitrary transformations of 
the parameter r then, effectively, it selects a specific r from all the possibilities and defines it 
as physical time - that is, as that definition of time which is consistent with our constraint on 
disc forces. 

Thus, at face value, we now have two equations for the system - one from the variational prin- 
ciple (not shown) and (jSj, above. However, it transpires that if one assumes that light traces 
mass in a disc galaxy (that is, no significant dark matter) then, typically for light-inferred mass 
distributions, one finds 

which corresponds to density distributions which fall off more quickly than 1/R but no 
quicker than 1/R^. Since B, defined at (|2||, is necessarily positive for any M, then equation ^ 
is necessarily complex which, in turn, implies K = q + pj for real parameters {p, q). Thus, in 
the final analysis we end up with three equations - which implies that the theory determines 
the dynamics and the mass distribution for a disc galaxy. 



4.2 The power-law as a limiting velocity distribution 

The same theory that gave almost perfect modelling for the sample of eight LSBs considered 
in Roscoe [2J (astro-ph/ 0306228 ) (data provided by McGaugh^) also provides an understand- 
ing of why the simple power-law model provides such a high-fidelity statistical resolution of 
middle-disc ORG data over large data sets. 

In effect, we find that the power-law solution, V = AR", is the only solution which is bounded 
over finite discs and that all other possible solutions become rapidly unbounded and there- 
fore, by definition, cannot persist over substantial radial ranges within stable discs. However, 
according to the theory, these non-power-law solutions can persist over small radial ranges 
and are therefore, in terms of the theory, responsible for much of the complexity commonly 
observed within optical discs. However, since, by definition, stable discs have bounded ve- 
locity distributions, then it is the power-law which is the dominant underlying form. 
To obtain a quantitative understanding of these remarks, we must consider the basic theory: 
the density-distribution equation derived in Roscoe \21 (astro-ph/ 0306228 ) can be written as 



(6) 



where e = -1 f- and = 1 -R% p = — 

p dR mi 2'kR oR 

where p is surface density, A4. is the basic mass function of the theory, mi is a scaling parameter 
having dimensions of mass and (j9, q) are free dimensionless parameters. The corresponding 
rotational velocity equation can be written as 



I dV _l 
VdR ~ R 



1 , 2\ P ( p^ + q \ 



(7) 



^Department of Astronomy, University of Maryland, USA 
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It is easily shown that if 

M =mo\n doR (8) 

where tuq has dimensions of mass and is dimensionless, then p ~ 1/ R^, = 1 and \I/ = 
const etc with the final consequence that the velocity equation gives V = AR" exactly But 
also, by the considerations of ^4.11 we know that p cannot fall of more rapidly than 1/R^ - 
otherwise equation ^ becomes real, and the whole analysis becomes invalidated. We are 
therefore left to consider solutions for which p~l/i?^, k < 2: from ^ it is easy to see that, 
for such solutions, ^ as i? increases. But, by Q, dV/dR ^ oo as ^ ^ so that, as stated, 
all other possible solutions are unbounded over any finite disc. 

To conclude, the power-law solution is the only solution which is bounded over a finite disc 
and therefore, according to the theory, provides the underlying dynamical 'shape' of any 
stable disc. 



5 The Development of the quasi-relativistic disc-model 

We show, in appendix 13 that it is impossible to have non- trivial power-law solutions merg- 
ing into flat solutions within the context of the quasi-classical one-clock theory. The solution 
to this problem turns out to lie in the quasi-relativistic generalization of the quasi-classical 
model: Specifically, it transpires that the flat component of rotation curves is defined in a hy- 
perbolic space which arises directly from the quasi-relativistic generalization whilst the opti- 
cal mid-region component remains associated with the quasi-classical model. The boundary 
between the two components can be considered analogous to a shock front, with the transi- 
tion across this boundary being managed by appropriately defined jump conditions. 



5.1 The Euler-Lagrange component 

The basic definition for the metric tensor, given at (H)) remains unchanged - except that now, 
the indices vary over 1..4 where = ct and t represents conventional clock time and c is 
a scaling parameter with the dimensions of velocity. In cylindrical polars, the line element 
becomes 

ds^ = gadx'dx^ = -^^dr^ + — — r^dO^ + 2——drdt + ^^dt^ 
or^ \r or J orot ot^ 

which, in dimensionless form, is given by 

where R = t/tq and T = t/to where (rg, to) are scaling constants, as before. Now introduce 

9i?2 ' RdR' dRdT' ' 

noting that (A, B, C, D) all have dimensions of mass, so that 

ds"^ = AdR^ + B R^de^ + DdRdT + EdT"^ (10) 

from which we get the variational principle 

1/2 



1 = j Cdr, C = ./g~¥i^ = (aR^ + BR^d"^ + DRT + ET^) 



(11) 



9 



As before, this is homogeneous degree zero in the parameter r, so that the equations of motion 
form an incomplete set given by 

^ {^^^Z } = ^ ^^^^ = (^^^ + BR^e^ + DRT + ET^^ ^'"^ (12) 

where mi is a scaling constant with dimensions of mass, together with 

In a classical context, the equation corresponding to (|T2)l would give angular momentum 
conservation. Clearly, angular momentum is not conserved here (the forces are not assumed 
central) - but something is. For convenience, we shall refer to ((T2)l as the equation for the 
conservation of generalized angular momentum. 



5.2 Physical time and mathematical closure 

As before, the problem is that the system is invariant wrt to arbitrary (monotonic) transforma- 
tions of the temporal ordering parameter, r, which implies that physical time is not defined 
- some physical condition is missing from the system. To progress, we could adopt the GR 
solution of introducing the notion of particle proper time - but that would entail ignoring the 
considerable success of the earlier approach (above and Roscoe [2], astro-ph/ 0306228 ) which 
depended critically upon the idea that the ratio of the magnitudes of transverse to radial 
force components was constant within any given disc. We therefore seek to generalize the 
approach of §4. II to the present case using the line element ((TO)l . In the following, we sketch 
out the analysis, giving the details in appendix |Bl The only factorization of the line element 
(ITOb similar to that of (Hb turns out to be 



ds^ = ds-ds where ds = SiVAdRR + S2VB Rd9 9 + S^Ve dTK, (14) 
where R and 9 are unit vectors as before, subject to the condition on Ai that 



D = 2SA-/AE^{^—j (15) 

We can now use ((T4)l to define a generalized velocity, ds/dr = V, and a corresponding gen- 
eralized acceleration, dV /dr = V, in the space in a way which is directly analogous to the 
development of ^4.11 The arbitrariness associated with the temporal ordering parameter, r, is 
removed by applying the physical condition that the ratio of radial to transverse generalized 
accelerations is to be constant within the disc. Adopting the conventions that 



then this physical condition gives rise to the additional equations 

^{^Vg) = -q^}l-p^f^SoV^Vn + SrV^f} (16) 

^(5ov^y« + 5iv^T) = -p./B^ + q^{So^/^VR + S^^/^f] (17) 
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whilst (ITSb becomes 

In the above, Vg = R9 and Vr = R are the rotational and radial components of the generalized 
velocity, V, and (p, q) are dimensionless parameters directly analogous to those of §4.11 
To summarize, our equations of motion are given by (|T2)l. ((T3)l , (O and ((T7)l together with the 
condition dTSt which must be also satisfied. Thus, we now have five equations to be solved 
for four unknowns, {R, R9, T, M) so that the system appears to be overdetermined - unless 
there is redundancy somewhere. We shall show that this is, in fact, the case. 

5.3 Linear dependencies with the system 

We show, in appendix O how (|T2)l. (O and (|T7)l lead to the equations 

y/B^ve - SoV^vr = SiV^ (19) 

Rd^ i?(9^ 
where 



ve = -r, vr = -v-, = S2\ I , ^2 = ±1 (21) 

T T y mi 

The structure of ((19)) and ((20)l is fundamental to the system. By removing the time-dependent 
terms, it is easily seen that the quasi-classical form of this latter pair is given by 



-}ve~-r.^VR = 



from which it is immediate that non-trivial solutions for vg and vr can only exist if 

This equation is, in fact, a rearranged form of the quasi-classical mass-equation, ©, which 
immediately suggests the possibility that (IT9|) and (|20)| are similarly linearly dependent. If so, 
we must necessarily have the conditions 



(1 + ^) 

q"^ 



simultaneously satisfied. We immediately see that this is only possible if the additional mass 
condition 

S.V^_ = S„V^^ (23) 

holds - remember that A, B, D, E, ^ are all defined in terms of j\4. But j\4 must also satisfy 
(IT8|) , so we must now consider whether or not such a thing is possible. 
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5.4 The Mass Equation 

Equation ((T8)l has two basic classes of solution, given by a logarithmic class 

M{R,T) = moHl + R + T) + do, (24) 

where (mo, do) are constants, and a variables separable class. The constant do is irrelevant here 
since only derivatives of Ai appear in the theory; we ignore it from now on. The variables 
separable class, which is considered in appendix O for completeness, turns out to be associ- 
ated with degeneracy in the system and so plays no part in the main body of this paper. The 
logarithmic class concerns us directly and we ask under what conditions (if any) the solu- 
tion (f24)l can satisfy condition (|23)l . Using the definitions © and (|2T|) we easily find that the 
logarithmic class can only satisfy (|23)l if the condition 

is satisfied. Noting that R = —{1 + T) is excluded by the structure of (|24|l. then we must have 

R=l + T, SoSi = -1. 

With this latter condition, we find that M{R, T) = nio ln(l + i? + T) is a solution subject to 
R = 1 + T and the conditions (f22)l which now reduce to the single algebraic equation 

Mq-S,2'''\p{2-^^+S,{l+p' + q^)sf^^ = (25) 

where M = mo/ mi, Si = ±1 and 5*2 = ±1. To summarize, the solution can be expressed as 

M = mo\n2R (26) 

on the characteristic R = 1+T subject to the condition that the dimensionless mass parameter 
M satisfies (|25||. 



6 Rotational velocity = constant in the hyperbolic disc 

In this section, we obtain the extremely interesting result that, when the mass distribution is 
as described in §5.41 then the rotational velocity is constant with a value determined by the 
parameters {p, q) and the algebraic equation (|25)l . The four basic (but rearranged) equations 
of motion are: 

yEWe-Sov^VR = SiV^f (27) 

A{_L(,,,...)}_| (.s) 



(I+P^Kp+^I Rd^ Rd^. 
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The first two equations arise directly from the variational principle ((TT|l . whilst the third and 
fourth equations arise from the conditions ((16)) and ((T7)l . In ^5.31 we showed the subsidiary 
mass condition ((T8)l could be satisfied by making ((27)l and ((30)l linearly dependent so that, 
in effect, we can consider the latter of these to be redundant. Consequently, we have three 
equations ((27|) . (l28b and (l29b say, for the determination of {Vq, Vr, T). 

To progress, (l28t must be expanded explicitly and this is done in appendix |El where we show 
how it reduces to an algebraic equation in the dimensionless velocity vg with the solution 

ve = . ^ (31) 
V2-M 

where M = mo/mi must satisfy dZSb . That is, the only possibility for rotation velocities in the 
disc of the quasi-relativistic theory is that of constancy - the rotation curve is flat. 



7 Jump Conditions 

We are concerned with the transition from the quasi-classical rising part of the rotation curve 
to the flat quasi-relativistic disc and, in all that follows, we use the superfix * to denote quanti- 
ties in the quasi-relativistic disc. It transpires that everything reduces to deriving the relation- 
ships between the quasi-classical parameters (mo, mi) and their quasi-relativistic counterparts 
(mo, mX) and, to this end, we use jump conditions which quantify 

• conservation of generalized angular momentum - cf equation ((T2)l : 

• conservation of mass-flow across the transition boundary; 

• no-slip for the rotational flow at the transition boundary. 

The first of these requires some discussion: within the context of the present model, forces 
within discs are not central forces - consequently, angular momentum is not conserved. How- 
ever, if we consider ((T2)l we see that the quantity B?6B / C (which plays the part here of clas- 
sical angular momentum) is conserved in both the quasi-classical and hyperbolic discs. Ac- 
cordingly, it is this quantity which must be conserved across the transition boundary. 



7.1 Conservation of Generalized Angular Momentum 

Reference to ((T2)l shows that the quantity Bj^OB / C must be conserved across the radial tran- 
sition boundary, R = Rt say. Using the no-slip condition Vg = vq at R = Rt then this 
requirement reduces to 

B* _B 

across R = Rt ■ Hence, following appendix|Fl we obtain the condition for generalized angular 
momentum across the transition boundary as ml = mi. 



7.2 Conservation of mass flow 

The mass-flow conservation across the transition boundary, R = Rt, can be stated as: 

P*'"*r\r^r^ = PVr\r=r^ ■ 
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Hence, following appendix O we obtain the condition for mass conservation across the tran- 
sition boundary as 



<o*c V2-M* _m* 



^ = ^S*,So^ {S;-Sl), M* = ^ (32) 

mo 



^ = S2VT^, M = — (33) 

mi 



together with 



At this stage, we can note that ^ 7^ necessarily otherwise, as is easily seen by reference to (l39t 
and ()40)l . dynamics in the quasi-classical disc becomes indeterminate. It follows immediately, 
from above that 

S^s; = -1. (34) 

7.3 Summary of jump relations 

For convenience, we collect the jump conditions in the form we use them, below: Using (|34|l 
then (|32)l and (|33)l give, after squaring. 



^ f2-ivr) = i-r^ 



2 

mo / ■ ■ V"^i 

where we remember that M* = ml/m\. Use of the generalized angular momentum jump 
condition m\ = mi now gives, after a bit of algebra, 

M'm*- f^V + 2-M*=0 (35) 

which we refer to as our first jump condition. Use of (l34t in (l32b gives 

^ = ^^*^o^iV2 - M*, M* = ^ (36) 

mo mi 

which we refer to as our second jump condition whilst our third is simply 

mi = mi. (37) 



8 Forbidden zones in the parameter space 



In the following, where the details of merging the rising part of a rotation curve (typically, 
the optical rotation curve) with its flat radio extension are considered, we remember that {p, q) 
are dimensionless disposable parameters which appeared firstly after ^ via K = q + pj and 
similarly in appendix |H1 We assume that they are fixed across the whole disc. 
To simplify matters, we represent the rising part of the rotation curve purely in terms of the 
power law solution V = AR". The theory says nothing about A since this would generally 
be determined by initial conditions, but is explicit about the structure of a in terms of the free 
parameters in the theory. Specifically, by ^ and the sentence which follows it, we have 



a 



p fp^ + q^ 
q \ 1 ' ^ 



(38) 
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and we shall consider the constraints placed upon this quantity by the theory Note that, in 
practice, we are only concerned with < a < 1 since there are essentially no optical rotation 
curves outside of this range (Roscoe | |27ll30t ). 

The constraining equations are, from ^ in the case of a power-law solution on the quasi- 
classical disc, the condition 

together with 

{/ M*\ I M* 1 

p[2- — ]+s;{i + p' + q')^i- — \=o- sis; = -i, 

from ()25)l on the quasi-relativistic disc. Together with (|35)l and (|36)l. these represent four con- 
straints on the five unknowns {p, q, M*, mo/rriQ, \E'). 

8.1 Brief computational details 

There are five unknowns and four constraints together with various unknown signatures. A 
simple (inefficient with lots of redundancy!) approach is as follows: 

• choose q on the range (—0.5 < q < +0.5) tna uniform random way; 

• choose S'J' = ±1 in a uniform random way, and then set 5*2 = —SI; 

• choose ^0 = ±1 and = ±1 in uniform random ways; 

• solve the four constraining equations for (p, M*, mo /wig, \E'); 

• calculate a from (l38t for each signature in turn and choose the value that satisfies < 
a < 1, if one exists, and record it; 

• repeat a large number of times. 

8.2 Results 

FigureHh plots the determinant of the Jacobian of the system against a. It transpires that there 
are four distinct branches of solutions, of which branch I and II are shown as solid lines and 
branch III as a dotted line. Branch IV is out of range of the diagram. The most immediate 
relevant features are: 

• on branch I, the Jacobian is singular at about a = 0.43; 

• on branch II, the Jacobian is singular at about a = 0.40. 

Thus, if the theory does describe aspects of reality for disc galaxies, then we would expect 
there to be a relative under-abundance of disc galaxies in that part of the parameter space 
which is in the neigbourhood of these singularities. 

However, there is a problem associated with considering a-data directly: whilst we know that 
a and In A are extremely strongly correlated (^ |27ll30J ) we also know that estimates of a from 
the data are an order of magnitude less accurate than those of In A (Roscoe |28 1). Consequently, 
rather than consider the a data directly, we consider it through the proxy of In A. Figure [it' 
shows what we find: 
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Figure 1: The left panel plots the determinant of the system Jacobian against the parameter 
a. It indicates the existance of singularites at a = 0.4 and at « = 0.43. The right panel, which 
primarily shows the In A distribution for Mathewson et al data, indicates the estimated posi- 
tions of these two singularities in In A space. They clearly occupy a region of the parameter 
space which is associated with an under-abundance of objects. 



• The figure primarily shows the strong four-peak structure in the In A distribution dis- 
cussed in detail in Roscoe [30J in the context of the discrete dynamical states phe- 
nomenology; 

• Low values of In A in the figure correspond to dim objects whilst high values of In A 
correspond to bright objects. Thus, there are virtually no objects in the ranges lnA< 3.5 
or In A > 5.5; 

• The mean a-values corresponding to each of the four major peaks are indicated in the 
diagram; 

• The vertical dotted lines of figure [T]^ indicate the estimated In A positions of the two 
singularities shown in a-space in figure [l^; 

• We see that, whilst they occupy the middle of the In A range where we might expect 
to see the highest density of objects, they in fact lie in a region where there is a very 
significant under-abundance of objects. 

We can conclude that the theory's predictions about the existence of singularities in the pa- 
rameter space is in strong concordance with that which is observed. 



9 Summary and conclusions 

The quasi-classical one-clock model has enjoyed a great deal of success within the context of 
dynamics in the optical disc - but it fails to admit the possibility of a transition to flatness 
within discs. We then considered the possibility that a generalization of the quasi-classical 
model into a quasi-relativistic form might provide a resolution to this problem and found this 
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to be the case. However, in order to merge solutions between the quasi-classical and quasi- 
relativistic parts of the disc, it then became necessary to consider jump-conditions relating 
the two components of the disc. It transpired that the equations representing these jump- 
conditions became degenerate in exactly that part of the parameter space within which disc 
galaxies were already know to be significantly under-abundant - that is, in the valley between 
the two peaks corresponding to In 4.17 and \nA^ 4.72 in figure [TJ 

Although the valley between 4.17 < In A < 4.72 is the deepest of the valleys in figure [H 
it is not unique - if the theory is able to provide an explanation for this deepest valley, we 
should also expect it to provide an explanation for the valleys between 3.90 < \nA < 4.17 
and 4.72 < InA < 5.10 - either in explicit terms of why these regions are underpopulated in 
discs or in explicit terms of why discs appear to prefer the peaks. As things stand, the two- 
component disc model does not do this. However, there are two major ways in which this 
latter model is significantly over-simplified: 

• Real disc galaxies have a central bulge region. In modelling terms, this could be rep- 
resented by an interior rotating spherical distribution so that the over-all model would 
be a three-component disc model which would have two sets of jump-conditions placed 
upon the dynamics - one from interior sphere to the middle disc and one from the mid- 
dle disc to the exterior flat disc. We might then expect these extra jump-conditions to 
provide additional explanations for observed features; 

• Our model assumes perfect rotational symmetry within discs whereas real discs have a 
very strong tendancy for a 180° symmetry - that is, to have two symmetrically opposed 
spiral arms. We might expect the additional complexity arising from accounting for this 
two-arm structure to give rise to significant improvement in the explanatory power of 
the theory. 

To conclude, the theory, which has the underlying structure of a LeSage-type vacuum gravi- 
tation theory, has 'clocked up' sufficiently many significant successes that one is inclined to 
believe that, at some level, it is reflecting reality. This alone is enough to justify a continued 
effort in its development. 

A The possibility of transition to flatness within the quasi- 
classical model 

It is sufficient to consider only equations ^ and ^ for which the parameters (p, q) are in- 
dependent and, once chosen, fix the space of possible solutions. Suppose we consider the 
general power-law solution which necessarily arises whenever p = k/ R^. In this case, © 
gives 

e2 = l, p^i^ + + q^)^i +p = Q (39) 

so that, for any choice of (p, q) then ^ is defined to have one of two possible values, both 
constant. Thus, ^ becomes the simple power-law equation 

I dV _a 
V~dR ~ R 

for some constant a = a{p, q) which gives the exponent in the solution, V = AR"" and, in 
Roscoe 1 30 1, we have seen that it is this class of solutions which defines the mean behaviour 
of optical rotation curves. 
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Now, clearly, a flat rotation curve is a special case of the power-law for which a = 0. Thus, in 
this special case we have the additional relation 



„.l^£^(£!±^) 1 (40) 

q \ q J ^ 

where the =f arises because = 1. Thus, (|39)l and ()40)l together define a relationship /(p, g) = 
between the parameters (p, q) which must hold if a = 0. It follows that if the conditions for 
transition to flatness in the radio disc are satisfied within any given disc (that is, /(p, q) = 0), 
then the conditions for a non-zero power-law solution within the optical disc are excluded. 
But this is contrary to the conclusions the analyzes of Roscoe [27, 30] for the generality of 
optical discs - that is, there are substantial reasons for believing that the transition to flatness 
does not occur within the context of the quasi-classical model. 

B Physical time and mathematical closure 

The omitted details for the analysis of ^5.21 are given below. Using the line element ((TO)l we 
find 

ds^ = ds-ds, ds = Si\^dRR + S2VBRdee + S3VEdTR (41) 
subject to the condition that M must satisfy 



We can now use (HTJ to define a generalized velocity in the space 

V= (^Siy/AR + S3^/Ef^R + S2VB Ree 
from which a generalized acceleration vector, V, can be defined as 

V = |5i^(v^i?)+53^(v^T)-52V^i?^2}R 



(42) 



+ !^S2^(yBRe^ + SiVaR9 + SsVEfe"^ 9. 

Using the condition that the ratio of radial to transverse accelerations is to be constant within 
the disc, then, after multiplying through by S2, this latter equation leads to 



i^S,S2^(yAR)+S2Ss^[VEt)-VBRe'^ = 

K 1^ (VbRO^ + S1S2VAR9 + S2S3^fe^ (43) 

as the generalization of Since, from the earlier analysis, we have that A < 0, it follows 
from (|42)) that E < also. Since B > necessarily, then ()43)l is also complex so that we must 
assume that K is complex, as before. Adopting the conventions that 
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and defining 

Ve = Re, Vr = R, K = q + pj 

then (|43)l gives 



R 

whilst (I42I1 becomes 



D = -2S.S,V^V^ ^ =W^- 

C A Rearrangement of the equations of motion 

The key to a successful resolution of the system lies in recognizing redundancy in a certain 
subset of the four basic equations of motion. This subset is obtained by processing equations 
((T2)). (IT6t and ((TTjl in a certain way: From ((T2)l with use of (ITSt we get directly 

V^B{BR^ - mi) = mi {aV^ + DVrT + Ef^^ 

= mi (^SojV^VR + SijV^t 
i 

Defining 



^ = ^2^/1 , 52 = ±1 

mi 



then we get 

Vg/B^ = SoV^Vr + SiV^f. (44) 
We can now write equations (O and (|T7)l as 

We now eliminate d{VBVg) / dr between these equations to get, after some work. 
Finally, using the identity 

^ = — = — -I T 

dr dR dT 
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(there is no 9 dependency) we find 



^+ — — r-^m^^=^&f^ 

Noting that (HH) and this latter equation are homogeneous algebraic relations between Vg, Vr 
and T and that the ordinary transverse and radial velocity components are given vg = Vg/T 
and vr = Vr/T, then we can write them as 

Rd-^ Rd^ 
The structure of these two equations is fundamental to the system. 

D The Variables Separable Class 

The variables-separable class solution of equation (ITSb has the general structure 

M{R, T) = mo{R + PoY{T + A)^-'^ + /32 
where (mo, (3q, (32) are constants. 

We ask under what conditions (if any) this solution (with /5o = /?i = for simplicity) can 
satisfy condition (|23)l . Using the definitions Q and (jZU we easily find that the variables 
separable class can only be a solution for the particular case M{R, T) = tuq^/ RT, and Sq = 
Si = ±1. If these conditions are satisfied, then we find that both of ((22)l reduce to a single 
algebraic condition having the general structure 

f(RT, M) = 0, where M = — 

mi 

and mi is the mass-scaling constant in the function \E' defined at (f2T|l . But f{RT,M) = 
implies directly that RT = const so that, finally, we find that the variables separable class is a 
solution subject to the conditions 

mo 



M(R,T) = VrT, RT = const, SqSi = +1, /(I, M) = where M = 

mi 

This corresponds to = const on the characteristic surface R = const /T. As we shall show, 
it transpires that this class of solution is associated with a degeneracy in the system which 
removes them from the current discussion. 



E Expansion of an Equation of Motion 

Equation (|28)| is given by 

iL|_L(,,,...)}_|.o (45, 
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where 

L = (^AR^ + BR^e'^ + DRT + ET^^ ^'"^ = f {Avj, + Bvl + Dv^^ + E) ^''^ = TL 

and A, B, D and E are defined at Thus (|45)l can be written as 

d f ^ ^r. ^xl 1 fdA 2 dB ^ dD dE\ 

dT 1 2L + - 2L (^^^ + ar^^ + ar^^ + j ^^^^ 

where we have used the identity 

1 d _ d 
fdr ^ dT' 
Equations ((T2)l and (IT9b give, respectively 

L = -^=RBvg, vr = So ^{R^ ^ve - S'ogi ^'C^ . 

Substituting these in (|46|) , we find that the velocity derivatives, which all arise from the term, 

d [Dvr\ 
dT\ 2L J 

vanish identically so that (|46)l becomes purely algebraic in the velocity terms: in terms of L, vr 
and vg it becomes 



M 2 dB . ^ d ^/B \ dE dE 



using the definitions of A, B, D and E given at and using Maple, we obtain the remark- 
able result that 

1 

ve = ± , = 
V2-M 

on the characteristic R = 1 + T where M = mo/mi. 

F Conservation of Generalized Angular Momentum 

In the following, we give the details omitted in §7.11 Reference to ((T2)l shows that the quantity 
R^OB/C must be conserved across the transition boundary. Using the no-slip condition Vg = 
ve at R = Rt, this requirement reduces to 

B* _ B 

at the boundary so that, using ((TT|l for C and the corresponding definition for C* we get 

B* _ B 

\/A*v*^ + B*vf \/Av\ + Bvj + Dvr + E ' 
Using (jTSjl, (llSb and (Il9t this latter condition becomes 



^vf - v*ji ^-v\ + 2vl - 2SoSiVR - 1 
We now consider the left-hand side and right-hand sides separately: 
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Left-hand side 

Using (150) 1 and (l33) l we get 



Ihs = , ^ ° = ^ ° = ^ = = v/^V2-M 



where the no-slip condition, f g = ve, and dSTl have been used for the last step. 
Right-hand side 

Equations (ISTJ and (|3T|l together give the scaled radial velocity in the hyperbolic disc as: 

Vr = So{S2 - Si). 

Using (l34t , we deduce immediately that vji = —2SoSi. Using this and dSTl) we get 



rhs = — = = \/rrin — = x/miv 2 — M. 



lhs=rhs 

Putting these together, we finally get 



mi = 777-1 



as the jump-condition which guarantees the conservation of generalized angular momentum 
across the transition boundary. 



G Conservation of mass flow 

In the following, we give the details omitted in §7.21 If we use R = Rt to denote the radial 
transition boundary, the mass-flow conservation across this boundary can be stated as: 

P*'"*r\r=Rt = P^r\b=Bt ■ 

Since 

on both parts of the disc, then we can write the mass-flow condition as 

B*v*R\R=n, = Bvn\n=nr- (47) 

Thus, we need appropriate expressions for (B*, t>|j) and (B, vj^). From §4.21 we have already 
noted that M* = mg In i? so that using the definitions given at then 

A* = B* = ^ (48) 

whilst in §5.41 we have noted that M = moln(l + R + T) on R = 1 + T so that using the 
definitions given at ©/ then 

mo r._ ^0 p _ mo 
^ - -4^' ^ - 2^' ^ - ^^^^ 
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also on R = 1 + T. Noting that the quasi-classical relations can be obtained from the quasi- 
relativistic relations simply by dropping of all time-dependent terms then (|27)l gives, after 
using dlHll and (EUl, 

v-k = s;^^^*v; = s*^*v; (so) 

whilst use of (|27|) directly gives, after using (|49)l and (|2T|) . 



vr = So — - ^ve - SqSi — = S0S2V2 - Mvg - SqSi. (51) 

Hence, with these definitions of B, B*, vr, v*p, and the no-slip condition vq = Vq at transition, 
then (|47)l can be written as 

'^S*^*ve = ^ (^SoS2V2^1dvg - SoSi 
i 



If we now use the flatness solution given at (|3T|I we obtain our jump-condition for mass con- 
servation as 

ml 2 mi 

where also, given the definition at (|2T|) . 

mn 



^* = S*Vl-M*, M* = 



ml 
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